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We explore the finite temperature phase diagram of the anisotropic XY spin chain using the 
Quantum Chernoff Bound metric on thermal states. The analysis of the metric elements allows to 
easily identify, in terms of different scaling with temperature, quasi-classical and quantum-critical 
regions. These results extend recent ones obtained using the Bures metric and show that different 
information-theoretic notions of distance can carry the same sophisticated information about the 
phase diagram of an interacting many-body system featuring quantum-critical points. 
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I. INTRODUCTION 



Quantum Phase Transitions (QPTs) are characterized 
by a dramatic change in the ground state of a many- 
body system prompted by a change of the parameters 
{A} specifying the hamiltonian H({\}) of the system. 
Unlike classical phase transitions, which are driven by 
thermal fluctuations, quantum phase transitions occur 
at T = and are driven solely by quantum fluctuations. 
In particular, second order quantum phase transitions 
are characterized by a gap between the ground state and 
first excited state that vanishes at criticality [1J. 

The analysis of these transitions has benefited from 
tools of Quantum Information theory. The von Neumann 
entropy and fidelity applied to many-body systems can 
identify phase transitions and reveal different scaling be- 
haviors at different regions of the phase diagram |2]-[7]. 
More recently it has been shown that the quantum fi- 
delity - a distinguishability measure between quantum 
states - can identify the quantum phase transition by 
comparing two ground states corresponding to slightly 
different values of the coupling constants {A} [8]-|22j. 
This new approach provides an alternative to the study of 
phase transitions using order parameters and symmetry 
breaking patterns, which depends on a priori knowledge 
of the physics of the problem [23] . 

Another distinguishability measure for density opera- 
tors was very recently introduced: the Quantum Chernoff 
Bound [24] [25] . Suppose there are n copies each of two 
density matrices p and a, i.e., p® n and c®", and one set 
has been given to us. We know what p and a are, but 
we don't know which set we received. There is always a 
two-outcome POVM {Eq,E{\ that measures and distin- 
guishes one set of n-copy states from the other such that 
the probability of making a misidentification, P e , is min- 
imized. In general, this probability will depend on the 
number of copies. In the asymptotic limit of n very large, 
the probability of misidentification has a particularly 
simple exponential dependence on n given by P e , m in = 
e -«?QCB _ Thg quantity £,qcb is called the Quantum 



Chernoff Bound. It is a function of the density matrices 
only, with £,qcb = -log (min < s <i TV(p a <7 1-8 )). The 
classical version of this problem was analyzed for the first 
time by H. Chernoff [26] . It took 50 years to prove the 
quantum version of this problem. 

The Quantum Chernoff Bound has many interesting 
properties. Among these is monotonicity under CP 
maps, which makes it a valid distinguishability measure 
[24 . In addition, £,qcb h&s an operational meaning aris- 
ing from a statistical inference problem. Intuitively, we 
see that for a fixed probability of error P e , the larger 
£,qcb is the smaller the number of copies of p and a we 
will need in order to distinguish them. 

In the present work, we apply the Quantum Chernoff 
Bound as a distinguishability measure to the manifold of 
Gibbs or thermal states associated with the anisotropic 
XY model in a transverse magnetic field. In the spirit of 
the information-theoretic and differential geometric ap- 
proach advocated in [10], we will compare states corre- 
sponding to hamiltonians with slightly different values of 
the parameters {A} to derive a metric tensor for the pa- 
rameter manifold itself. This metric detects the second 
order quantum phase transitions of the XY model and 
shows the influence of finite temperature effects over the 
zero-temperature critical points. It is shown that the dif- 
ferent regions of the {T, A} phase diagram of the system 
above the T=0 critical points can be characterized by the 
different scaling behavior of the metric tensor with tem- 
perature. This paper parallels and extends the analysis 
reported in [2T] for the Bures metric. While the Quan- 
tum Chernoff Bound is associated with quantum state 
discrimination, the Bures metric is related with another 
natural probabilistic protocol: quantum estimation [22] . 
Our findings about the Quantum Chernoff Bound metric 
show that the same information about the phase diagram 
of an interacting many-body system can be obtained by 
using two independent distinguishability measures. 

This paper is organized as follows. In Section II we in- 
troduce the many-body system we analyzed and obtain 
the corresponding metric tensor in Section III. In Section 
IV we carry out the thermal analysis of the metric ele- 
ments. A global property of the metric is discussed in 
Section V, while conclusions and further research direc- 
tions are outlined in Section VI. 
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II. THE XY MODEL 

We analyze the quantum phase transitions of the one- 
dimensional spin chain XY model in a transverse mag- 
netic field given by the hamiltonian 



where p = X^P«I*)(*I ^ s a spectral decomposition of the 
density operator. Since we are interested in the finite 
temperature XY model, we consider thermal states of 
the form p = e z . Taking the spectral decomposition 
of p over cigcnstatcs of the hamiltonian, we can split the 
metric into two parts so that ds 2 = ds 2 + ds 2 ncl where 



H = - 



E 



1 + 7 
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'i+l 



'i+ 



i + Aof, (1) 



where the total number of spins N is odd, A e S 
is the transverse magnetic field along the z-axis, and 
7 G [—1,1] is the anisotropy parameter. For 7 = ±1 
and 7 = we obtain the Ising and XX models, re- 
spectively. This hamiltonian can be diagonalized using 
Jordan- Wigner, Fourier and Bogoliubov transformations 
27J. The diagonalized hamiltonian is given by H = 



E 



k=i 



_! Afc6l6fc, with Afe the quasi-particle energies 



given by A fe = ^e 2 + A 2 , with e k = cos ( 2 j^) - A and 
Afc = 7sin(^^). One-particle excitations are created 

by the fermionic operators b k = cos(^ L )<ij c + isin(^ L )d_fe 
acting on the ground state 



MA, 7)) 



fe=i 



;cos(^)|00) fe ,_ fe 



•isin(y)|ll) fci _ fc ), 

(2) 

with dk\00)k,-k = d-k\00)k,-k = h\9s(\,i)) = 0, and 
cos(6k) = efc/Afc. This model exhibits a quantum phase 
transition at two regions in the parameter space {7, A}: 
at critical lines A = ±1, and 7 = for — 1 < A < +1. 
At those critical regions the system becomes gapless. In 
general the gap A is given by: A= |1 — |A||if|A| > |1— 7 2 | 

if |A| < |1 - 7 2 | (region 



(region A); A = | 7 |^/l 
B). 
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III. QUANTUM CHERNOFF BOUND 



dsi = 



ds„ r . = 



% 



(d Pl f 



tyj 



mdj)\ 2 ( Pl -p^ 



(5) 
(6) 



We call the first sum the classical part of the met- 
ric since it only depends on the Boltzmann weights of 
the density operator. We label the second sum the non- 
classical or quantum part, because it explicitly depends 
on the states. 

In passing we would like to notice that ds 2 nc does 
not, strictly speaking, define a metric in the parameter 
space. Indeed one can have curves t — » p(t) of den- 
sity matrices where all the relevant \j) are fixed, e.g. 
p(t) = ^2jPi{t)\i){i\- In other terms, the quadratic form 
defined at each point of the parameter manifold is just 
positive-semidefinite rather than positive definite and 
the metric matrix can have zero eigenvalues. This re- 
mark has to be kept in mind when one considers the 
zero-temperature limit (3 — > 00 where, as we will see, 
ds — > ds nc . Moreover, in this limit g nc is nothing but 
the Fubini-Study metric over the projective state space 
(the same is obtained by starting from the Bures met- 
ric). Once this latter is pulled back to the parameter 
space null eigenvectors can appear. These vectors corre- 
spond to directions along which changing the parameters 
results in indistinguishable states. 

To evaluate the first sum ([5|, note that 



The Quantum Chcrnoff Bound (QCB) is given by 

Zqcb = -log ( min Tr(/»V-)), (3) 

where P e ,min = e~ n ^ QCB is the minimum probability of 
error in distinguishing two preparations p® n and a® n in 
the limit n -> 00 [24] [25]. Defining Q{p,o) = e~^ CB , 
the quantity 1 — Q serves as a distinguishability measure 
between states. Notice that 1 — Q(p,p) = 0, and 1 — 
Q(p, a) = 1 if p and a are orthogonal. Considering two 
nearby states p and p+dp, this quantity induces a metric 
tensor on the manifold of density operators, with the line 
element given by 



ds 2 



1 



E 



\(i\dp\j)\ 2 



dpPi 



<Et-(E))- 



-PEi 



d^pi = -0(d y Ei-(d^E))- 
d xPi = -0(d x E i -(d x E))- 



z 

-0Ei 

~z~ 

,-pEi 



where E t = J2k n k^k, (n k G {0, 1}). 

Summing over states and incorporating the fermion 

statistics (n M ) = (1 + e /3A ")~ 1 and using the free-fermion 

(4) property (n^n^) - (n^)(n„) = S^ u (n^)(l - (n M )), we get 

the following six components for the 3x3 symmetric met- 



ric tensor defined by ds 2 = g^dx^dx", x^ S {(3, 7, A} 

" -E * 
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9%i 
9%x 

g 11 
<? 7 a 

9\x 



16^cosh(/3A fc ) + l 



A? 



1 



I67 ^ cosh(/3A fc ) + 1 



A^ 



1 



16 ^ cosh(/3A fc ) 



-efe 



/^ 2 



E 



I67 2 ^ cosh(/3A fc + 1) A 2 k 



E 



e.A 2 fe 



I67 ^ cosh(/3A fe ) + 1 A 2 



1 



16 ^ cosh(/3A fc ) 



lAf 



(7) 

(8) 

(9) 
(10) 

(11) 
(12) 



In order to find the metric g™J; corresponding to ds 2 nc 
we first use the expression for the eigenstates (pi to ob- 
tain 



d.s 



2 

nc 



zE ^*i:> g- 



4 ^ cosh(/3A fc ) + 1 



Differentiating 9^ along each of our three parameters 
gives the metric <?"£. This part of the metric has only 
three nonzero components, namely 






nc 
9~/X 



nc 
9 XX 



cosh(/3A fe )-l\ e 2A2 



4fMcosh(/?Ai) + i; Af 
cosh(/3A fc )-lNe fe Al 



1 y^ / cos 
h ^ Vcos 



47 



jE( 



h(/3A fe ) + i; A 



cosh(/3A fc )-l\ A 2 



cosh(/3A fc ) + l/ At, 



(14) 
(15) 
(16) 



Notice that g^a vanishes along the line 7 = 0. 



A. Quasi-Classical Region 

We now find the scaling behavior with temperature for 
the metric element <?",?, in the limit /3A >• 1. The scal- 
ing of all other metric elements can be obtained following 
procedures similar to those illustrated below. 

Since the system is gapped in this region, we have 
/3Afc > /3A 3> 1. This permits us to make the fol- 
lowing approximation in the integral expression for <?™5;: 
/ cosh(/3A fc )-i \ ^ ggAfc-g ^ i_ 4e -y9A» Then o 

approximated by 



™ can be 



s" 7 (T = 0) - 



2tt 



_ /3A /efesin(fc) 



A* 



The exponential e~ f3Ak behaves like a sharp peak in k 
centered around the absolute minima of Aj. and can be 
approximated further by expanding A^ to second order in 
k. The minima of A/, occur at k = for region A, A > 0; 
at k = ±arccos( 1 J x 2 ) for region B; and at k = ±7r for 
region A, A < 0. 

Due to the peak of e~^ Ak , the rest of the integrand can 
be approximated by its value at the minimum of Aj, in 
region B, or expanded until second order in k in region 
A. Taking the limits of integration to be from —00 to 00, 
this approximation results in Gaussian integrals, with the 
result given by 



'77 v 



(/?A » 1) = g%(T = 0) - f(Xn)T a e 



-A/T 



(17) 



where the exponent a of the temperature is equal to 3/2 
for region A, and 1/2 for region B. Each metric element 
will have a different /(A, 7). Similarly, the classical terms 
of the metric elements can be approximated by 



=(/3A»l)=/ l (A !7 )T° 



,-A/T 



(18) 



In TablelTlwe summarize the scaling in temperature for 
all the metric elements in the quasi-classical region. 



IV. THERMAL ANALYSIS OF THE METRIC 

ELEMENTS 

We now proceed to analyze the scaling behavior of the 
metric elements with temperature, for two characteristic 
regions in the parameter space {(3, 7, A}: quantum- 
critical and quasi-classical. For the quantum-critical 
case we analyze the scaling at the critical region. For 
the quasi-classical case we carry out the scaling analysis 
away from the critical region, for temperatures small 
enough so that the system looks effectively gapped, i.e., 
/3A ^> 1. We carry out the analysis in the thermody- 
namic limit, in which -^- — > k, ^~J fc — > =— J_ dk after 
rescaling g — > g/N. 



TABLE I: Temperature exponent a, for classical and nonclas- 
sical terms of the metric elements in the quasi-classical region. 
The nonclassical terms follow a behavior given by (171, while 
the classical terms behave like (18 1. 
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B. Quantum-Critical Region 

Consider now taking T —> at the values of 7 and A for 
which A = 0. Since in this limit (3 — > 00, all classical met- 
ric elements vanish due to the factor of (cosh(/3Afc) + l) _1 



in front. We are left only to analyze the nonclassical part 
of the metric, namely g™, g™^, and g^. There are three 
cases to consider: i) 7 = and A = ±1, ii) 7 = and 
-1 < A < 1, and hi) 7 ^ and A = ±1. 

Since g™^ and <7™£ are multiplied by an overall factor 
of 7, both vanish for cases i) and ii). 

The dispersion A& is an even function of k, so we can 
restrict ourselves to the interval k £ [0,7r]. For cases ii) 
and iii) A^ is linear in k about its root, whereas for the 
case i) A& is quadratic about its root. 

We now show how to calculate the scaling in temper- 
ature of the metric element <?™f; for case i). A similar 
procedure applies to g™^ and 5™^ for the other cases. 

The goal is to bound the metric element above and be- 
low by functions that have the same scaling behavior in 
j3. This will ensure that the metric itself must scale with 
the same exponent. 

Since the dispersion A/, is quadratic around the root 

,2 
k = 0, we can approximate it as A^ ~ ^-. Dehne the 

piecewise function 



f(P,k) 



/J 2 /,: 1 

16 
1 



for < k < 



for 



< k < IT, 



(19) 



which for all (3 and k satisfies 

f(f3,k) < cosh(^) 



1 



cosh(^) + l 



<f((3,k). 



(20) 



We can now split the integral into two parts 
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n cosh(/3A fe ) 
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At 



dk 
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cosh(/3A fe ) 



Ik 2 



1 e\ sin 2 k 



Air J_a_ cosh(/3A fc ) 

v7? 



At 



(21) 
(22) 



where we have taken A& ~ 4- and sin(fc) ~ k for the 
integral (21). Note that this is a good approximation 
for (3 — ► 00, since the upper integration limit becomes 
arbitrarily close to 0. 

We can now bound the hrst integral above and below, 



using the function (19 1 
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cosh(^) 



1 4 



C0S h(^) + lfc 2 



sTP 4 

dkfi^k)-. 



These bouri 



: ding integrals scale as /3 2 . Therefore the 
first integral (21 1 must also scale as /3 s for j3 — > 00. 



Over the entire interval [0,7r] we can bound A fe from 
above by -j^, upper bound sin (A) by k 2 , and replace the 



ratio of hyperbolic cosines in (22) by /(/?, k). This will 



bound integral (22 1 from above, as 



J^ r cosh(/3A fc )-le 2 sin 2 fc _£_ 

47r7^ cosh(/3A fc ) + l A* " " 4tt . _ 

~/3i 



f/A- 



3G 



Therefore, since the second integral ( 22 1 scales as a power 



of I 



or lower in (3, g™!r must scale as /Ja to highest order. 
Table [n] lists the scaling for all three non-classical metric 
elements. 






T -l 







-1/2 


const+0(T) 
const+0(T a ) 




T -l 



TABLE II: Scaling behavior with temperature at criticality, 
for the metric elements <?$a, g™1 and g"^. Note that g^\ and 
g"x are exactly for 7 = 0. 

A = ±l,7 = A = ±l,7^0 A£(-l,l),7 = 

nc 

9\x 



nc 

^77 



The behavior of the nonclassical metric elements in the 
quantum-critical region can be inferred from dimensional 
scaling analysis as was done in [21 , for zero temperature. 
Though carried out for the Bures metric, it applies also 
for the Quantum Chernoff Bound metric. From |21j we 
see that the scaling dimension of the nonclassical metric 
elements g%%, 3"^, and g™ is given by A„ c = A M + A„ - 
2z — d, where A M , A„ are the scaling dimensions of the 
operators that couple to /1, v £ {A, 7} and produce the 
quantum phase transitions, z is the dynamical exponent, 
and d the spatial dimensionality. The finite temperature 
in the quantum system transforms into a finite additional 
length dimension in the classical system so that 



9 



T 



A„ c / Z 



(23) 



The findings summarized in Table ITT] are consistent with 
this relation, assuming that: A 7 = 2, A^ > 3 (1st col- 
umn, z = 2); A 7 = 3, A\ = 1 (2nd column, z = 1); 



1, A\ > 2 (3rd column, 



1). The condition 



A x < d + z (d = 1) signals the relevance, in the renor- 
malization group sense [23] , of the operator weighted by 
the coupling constant x (x = 7, A). 

Intuitively, shifting some parameter towards a critical 
value should result in a significant change in the state 
of our system. That is, the direction in parameter space 
corresponding to maximal distinguishability should be 
towards or away from the critical regions. Indeed we see 
that the diverging metric element in the gapless quan- 
tum critical case corresponds to exactly this direction of 
maximal distinguishability. As an example, for A = ±1 
and 7^0 only g\\ diverges. Changing A and keeping 
7 constant here corresponds to the direction of maximal 
distinguishability, which is along the g\\ contribution to 
the line clement. In the Ising model limit, where 7=1, 
our results for the metric scaling agrees exactly with that 
of [2T], where the Bures metric was used. 



V. MAXIMUM EIGENVALUE 

In this brief section we consider the maximum eigen- 
value of the metric g for each point in parameter space. 
Since g is a symmetric matrix, this corresponds to the 
matrix norm \\g\\ = max|| v || =1 |(v,gv)|. The associated 
eigenvectors define the field of maximum distinguisha- 
bility directions. An analogous analysis for the Burcs 
metric has been presented in |21j . From the inequalities 
maxij \g%j\ < \\g\\ < 9maxij \gij\ one immediately sees 
that the maximal eigenvalue of g diverges iff at least one 
of the matrix elements of g diverges. In other words, ||<?|| 
encodes global information about the metric that one can 
extract without analyzing each of the matrix elements 
separately. 

Shown in Figure [T] is a contour plot of the maximum 
eigenvalues of g. The critical regions are clearly revealed 
as those values of the parameters for which some metric 
elements diverge for small temperature. 




FIG. 1: Contour plot of the maximum eigenvalue of g for 
-1.5 < A < 1.5, -1 < 7 < 1, and T = 1CT 2 . The color scale 
goes from (purple), to 3 (white). The critical regions can 
be clearly identified. 



Bound to study the finite-temperature phase diagram of 
the anisotropic XY spin chain. We have shown that the 
temperature scaling of the metric elements has a different 
form in the quasi-classical and quantum-critical cases. In 
the quasi-classical case, where the temperature T is much 
smaller than the energy gap A, the metric elements van- 
ish exponentially in 1/T; in the quantum-critical regime 
A/T <C l,T — > 0, the behavior is a power law in T. 
Moreover, in the quasi-classical case the specific behav- 
ior in temperature depends on the functional form of the 
gap A, and in the quantum critical case the nonclassical 
metric element which diverges in the limit of T — > is the 
one corresponding to the relevant parameter driving the 
QPT. We would like to notice that, in view of the oper- 
ational meaning of the Quantum Chernoff Bound metric 
in quantum state discrimination tasks, these results sug- 
gest that quantum criticality may provide a resource for 
quantum state discrimination protocols. This connection 
has been elaborated in [55] for the related quantum es- 
timation problem. Another problem that is worthy of 
investigation and has not been touched in this paper is 
the role of the curvature tensor associated to the met- 
ric g. The preliminary findings for the same quantity in 
the pure state case [TUj and in the Bures metric case 
[2T] suggest that the curvature might provide additional 
information about the nature of the critical lines and lo- 
cation of crossover regions. 

The analysis reported in this paper for the XY chain 
is quite close in spirit and technicalities to the one 
for the Bures metric in the Quantum Ising model case 
[5Tj . In both cases the important physical features 
of the low-temperature phase diagram of the systems, 
e.g., quantum-critical to quasi-classical crossovers, can be 
identified by resorting to a metric with a statistical distin- 
guishability meaning. This proves that the information- 
theoretic metrical approach to quantum criticality |10j 
is not bound to the use of a special single metric. Dif- 
ferent choices of distinguishability measure may provide 
essentially the same physical information. 
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